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Abstract 



^ . This paper is devoted to the study of time-dependent hyperbolic systems and the 

K^ I derivation of dispersive estimates for their solutions. It is based on a diagonali- 

sation of the full symbol within adapted symbol classes in order to extract the 
t:j- . essential information about representations of solutions. This is combined with a 

^N I multi-dimensional van der Corput lemma to derive dispersive estimates. 

in 

o, 

^D . Key words: hyperbolic systems, diagonalisation, dispersive estimates, oscillating 

time-dependent coefficients, van der Corput lemma 
1991 MSC: 35L05, 35L15 



1 Introduction 

Our aim is to consider Cauchy problems for hyperbolic systems 

DtU = A{t,D)U, f/(0,-) = f/o, (i; 
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where A{t, D) denotes a smoothly time-dependent matrix Fourier multipher 
with first order symbol 

A(t,0 gC°°(M+ xM^,C™'^™)nL°°(M+,5^(M")) (2) 

subject to certain (natural) assumptions which are described later on in detail. 
As usual we denote D( = —idt and we assume that eigenvalues of A are real 
modulo 'lower order terms'. 

The decay of solutions to the scalar higher order equations with constant 
coefficients of dissipative type was considered in [13]. Consequently, a compre- 
hensive analysis of the dispersive and Strichartz estimates of scalar hyperbolic 
equations of higher orders with constant coefficients, of general form, as well 
as for hyperbolic systems with constant coefficients, with applications to non- 
linear equations was carried out in [14] . In analogy to this, the results of this 
paper have the corresponding (by now standard) consequences for Strichartz 
estimates, and for the global in time well-posedness of nonlinear equations, 
which we omit here, and refer to [14] for the corresponding constructions. 

At the same time, equations of the second order with time dependent coef- 
ficients have been intensively studied in the literature, see e.g. [10], [11], [8], 
[19], [21]. In [6], the asymptotic integration method was developed for scalar 
hyperbolic equations with time dependent homogeneous symbols. The empha- 
sis there was placed on minimising the time-differentiability assumptions on 
the symbol, required by applications to the Kirchhoff equations ([5]). Indeed, 
the representation of solutions and dispersive estimates have been derived in 
[6] for equations of any order, by making assumptions on one time- derivative 
of the coefficients only, see Remark 21. However, the asymptotic integration 
produces an additional loss of regularity for high frequencies due to the fact 
that the amplitudes in the representation of solutions are symbols of type 
(0, 0). This is avoided in this paper since we construct the amplitudes as sym- 
bols of type (1,0). In [3] a special case of the results of this paper concerning 
differential hyperbolic systems with 7 = 2 (i.e., with non- vanishing Gaussian 
curvature of the Fresnel surfaces) has been treated. 

Our approach is based on diagonalising the (full) symbol of the operator, mod- 
ulo controllable terms, in order to gain a representation of solutions in terms 
of Fourier integrals. The basic idea of this diagonalisation scheme comes from 
the treatment of degenerate hyperbolic problems suggested in thee book of 
Yagdjian [22]. This representation is given in Theorem 20. Consequently, it is 
used to derive the dispersive estimates for solutions to the Cauchy problem 
(1) in Theorem 26 and Theorem 27, which are the main results of this paper. 
For this purpose, we derive more general estimates for oscillatory integrals in 
Propositions 24 and 25. The time decay rate of solutions to (1) is related to the 
geometric properties of the Fresnel surfaces of its characteristics in both con- 
vex and non-convex situations. In Section 4 we illustrate our assumptions for 



differential hyperbolic systems, for hyperbolic equations of the second order, 
and for general scalar hyperbolic equations of higher orders. 



2 Prerequisites 



The first step is based on a decomposition of the extented phase space M+ x 
M" into zones, a pseudo-differential zone containing bounded frequencies and 
a hyperbolic zone containing the large frequencies. It follows basically the 
treatment of [10], [19], going back to the papers of Reissig-Yagdjian for the 
treatment of wave equations with time-dependent propagation speed ([8], [7] 
and references cited therein, and in [10]) and the book of Yagdjian, [22], where 
the method of zones was introduced and applied to the investigation of the 
Cauchy problem for hyperbolic operators with multiple characteristics. 



2.1 Hyperbolic symbol classes 

We make use of the implicitly defined function t^ from 

(l + t5)|e| = iV(log(e + t^)r (2.1) 

with suitable constants A^ and v (usually required to belong to [0, 1]) and 
define the zones 

Zhyr>{N,v) = {(t,0|t > tj, Zp,{N) = {(t,O|0 < t < tj. (2.2) 

In Zhyp{N, u) we apply a diagonalisation procedure to the full symbol. The 
basic idea of this diagonalisation scheme comes from the treatment of degen- 
erate hyperbolic problems and is closely related to the approach of [10]. If 
u > 0, the hyperbolic zone is subdivided further into two parts, the regular 
and the oscillating sub- zones 

ZregiN, u) = Z^y^^N , 2u) Z,,,(iV) = Z^y^^N , v) \ Zhyp{N, 2u). (2.3) 

We denote the new boundary as t^. It satisfies (1 -|- t^)|^| = A^( log(e + t^) j . 
The parameter v is used to control the allowed amount of oscillations in the 
time-dependence of coefficients and amplitude functions. Based on the ter- 
minology from [22, Section 3.11.1] and of [7], [9], one speaks of very slow 
oscillations if z/ = 0, slow oscillations of z/ G (0,1), and fast oscillations if 
V = 1. The following definition of symbol classes and also the basic assump- 
tions introduced later on will depend on the parameter v. 



Definition 1 The time- dependent Fourier multiplier a{t,$^) belongs to the 
symbol class Sj^'j/{mi,m2} if it satisfies the symbolic estimates 



D^D"a(t,0 <C,,,|el 



1 / \u'^ m2+k 



1+t 



log(e + 1) 



(2.4) 



for all multi-indices a G N" with \a\ < ii and all natural numbers k < £2, 
where 

|ek = max(|e|,^(log(e + t))'). (2.5) 



If the symbolic estimates liold for all derivatives we write Si\i^t,{mi,m2} short 
for S^'^{mi,m2}. Furthermore, the definition extents immediately to matrix- 
valued Fourier multipliers. The rules of the corresponding symbolic calculus 
are simple consequences of Definition 1 together with (2.1)-(2.3) and are col- 
lected in the following proposition. 

Proposition 2 (1) Sj^'^'^{mi,m2} is a vector space. 

(2) S^i^f{mi -k,m2 + i} ^ S^j^'^'{mi,m2} for all £ > k > 0, £[ > £1 and 

(3) 5^';^' {mi, ma} ■ >S^'^'{m'i, m'2} ■-> 5^';^' {mi +m[,m2 + m'2}. 

(4) D,^D^5^;^^{mi,m2} ^ 5S;l"l''^^'{mi - \a\,m2 + A;}. 

(5) S'^li-1,2} ^ LfLliZregiN,u)). 



PROOF. The statements are simple consequences of the symbolic estimates. 
Since it is vitally important for us later on, we show the last one. It follows 
that 

I . .x|, / /•- (log(e + r))2- 

t Jt me + ry 



(log(e + r )) 



2u 



2v 



J t 



{\og{e + T)f --^ 

lel(e + r) 



-dr 



^ (log(e + t))^- ^ 



for any a G 5^y{ — 1, 2} uniform in (t, ,^) G Zreg{N, v). 



Of particular importance are the embedding relations (2). They constitute a 
symbolic hierarchy, which is used in the diagonalisation scheme, cf. Sect. 3.2.2. 
We define the residual symbol classes 



'HN,v{'m} = Pi SN,u{'m - k, k} 

k&Z 



(2.6) 



as intersections along the embeddings Sn,u{'>tT'1 — k, 1712 + k} ^-^ SN,u{'m'i, 1^2} 
for A; > 0. 

For later use we define the cut-off function Xpd(t,0 = x(l'Cl(l + t)(}og{e + 
t))'''/N) for X e <^o°°(^)' X(s) = 1/2 for s < f and x(s) = for s > 1 and 
similarly Xhypit,i) = 1 - Xpdit,0- 



2.2 Basic assumptions 

We collect our assumptions on the symbol A(t,,^). Throughout this paper we 
require 

{Al)f^f,^ Hyperbolic operator of first order with bounded coefficients. We 
assume that the matrix operator A{t, D) has a symbol satisfying 

A(t,Oe5^,fMi,0} (2.7) 

for some N > and u G [0, 1], which has all its eigenvalues in a strip 
spec A(t,0 ^Sc = {C eC:\lmC\<c} for all (t,0 e ZhypiN, v). 

We say the system is symmetric hyperbolic if we assume in addition 

ImA(t,0 = ^(A(t,0 - A*(t,0) G 5fc'H0,l}. (2.8) 

We write (Al^) for this stronger assumption. 

(A2) Uniform strict hyperbolicity up tot = 00. We assume that there exists a 
homogeneous of order one in ^ matrix Ai (t, ^) with A(t, ^) —Ai {t, OXhypiti ^ 
5^'^^{0, 1}, which has real and distinct eigenvalues 

Ai(t,o,...,A™(t,OG5^V'Hi,o}, 

the so-called characteristic roots of the symbol Ai{t.,^). In ascending order we 
denote them as Ai(t, .^), ..., Am(t, ,^). Furthermore, we assume that 

liminf min |Ai(t, w) - Aj(t, a;)| > (2.9) 

for all i ^ j. 

As a consequence of the first two assumptions we get 

Proposition 3 Assume (Al)^^^^^ and (A2). For all j = l,...,m the charac- 
teristic roots satisfy \j{t,^) G iSj^'^^{l,0} and for all i ^ j their difference 



satisfies (Aj(t, ^) — \j{t,^)) ^ G 5^'^^{— 1,0}. Furthermore, the eigenprojec- 
tion Pj(t,^) corresponding to \j(t,^) satisfies Pj(t,^) E Sj^'^'^{0,0}. 



PROOF. [Sketch of proof.] The properties of the characteristic roots follow 
from the spectral estimate \Xj{t,u)\ < \\Ai{t,uj)\\ together with the obvious 
symbol properties of the coefficients of the characteristic polynomial and the 
uniform strict hyperbolicity. The eigenprojections can be expressed in terms 
of the characteristic roots 

and again the symbolic calculus yields the desired result. 



For ai, . . . , Om G C we define diag(ai, . . . , a„i) to be the diagonal matrix with 
entries ai, . . . , a^ along the diagonal. For a matrix M, we define diagM as 
the diagonal matrix with Mn, . . . , Mmm at the diagonal. 

Proposition 4 Assume (Al)(,^^i^ and (A2). There exists a uniformly hounded 
and znvertzble matrix M{t, e S^j^;^^{0, 0} with M'^it, e >S^;^'{05 0} '^hich 
diagonalises the symbol Ai{t,S,), 

A,{t,i)M{t,uj) = M{t,0'D{t,0, V{t,0 = diag(Ai(t,0,...,A„(t,0). 

(2.11) 



PROOF. By Proposition 3 we see that the symmetriser of Ai{t,^) given 
as H{t,^) = EjPjit,OPji't,0 satisfies H{t,{) G S'^'f^'{0,0}. Our aim is 
to express M^^(t, ^) explicitly in terms of the eigenprojections Pj{t,^) G 
iS^'j,^{0, 0}. Let therefore be Vj{t,^) be (a smoothly chosen) unit vector from 
the one- dimensional j'-th eigenspace ran Pj(t,^) and M~^{t,C,)u be the vector 
with the entries given by the scalar products {vj, Pj(t, ^)u) = {vj, H(t, ^)u) for 
any vector -u G C". 

Note that Vj{t,C,) is unique up to a sign and locally in t and C, expressable 
as Vj(t,^) = Pj(t,^)e/\\Pj(t,^)e\\ for some fixed unit vector e chosen away 
from the orthogonal complement to the eigenspace. Differentiating this yields 
the desired bounds Vj{t, e >S^;^^'{0, 0} and therefore M-\t, e >S^;^^HO, 0}. 
The inverse matrix M{t, C,) has the vectors Vj{t, C,) as its columns and therefore 
M(t,0 e 5^;^'{0,0} as claimed. 

Remark 5 Under Assumption (Al^) the matrix M{t,^) can be chosen uni- 
tary, which simplifies some of the considerations later on. 



(A3) Generalised energy conservation property. We have to make a technical 
assumption which guarantees us that the lower order terms are not too strong. 
This can be written as 



sup 



lmFo{e,^)de 



< oo 



(2.12) 



for the matrix 



Fo{t,0 = dmg{M-\t,0{A{t,0 - Ai{t,0)M{t,0 - M-\t,0{^tM{t,0)). 

(2.13) 

Remark 6 This statement does not follow from symbol estimates as we only 
have Fq E S^^'J^^ {OjI}- ^^ ''^^^^ later use it in Theorem 15 relating it to 
uniform hounds on the micro-energy of solutions within ZregiN,!/) explain- 
ing the nature of (A3). Assumption (A3) is independent of the choice of the 
diagonaliser M{t,^), provided it satisfies the properties of Proposition 4- 

If we assume (A1+), then the diagonaliser M{t,C,) can be chosen unitary. This 
implies 

Proposition 7 Assume (Al^). Then the matrix M*(DtM) is self-adjoint and 
therefore 

Im diag M*(t, ODtM(t, = 0, (2.14) 

zn particular, Im Fo(t,0 = Im diag{M*{t,^){A{t,^) - Ai{t,^))M{t,^)). 



PROOF. Obviously = DJ = (DtM*)M + M*{DtM) and therefore 

(M*DiM)* = (DiM)*M = -{DtM*)M = M*{I)tM). 
But this implies that the diagonal entries of {DtM*)M are real. 



Assumption (A3) can be replaced by the weaker assumption that there exists 
a real scalar function S(t,^) G S{0, 1} such that ImFo(t, ^) — 6{t,C,)I satisfies 
the equivalent of (2.12), i.e. 



sup / ilmFoi9,0-Si9,Ol)de 

{s,0,{t,OeZreaiN,u) 

We refer to this weaker assumption as (A3^). 



< oo. 



(2.15) 



(A4) Weak dissipativity There exists a constant c and a function 7(t) such 
that 

ImA(t,0 + c|e|I+|7Wl>0 (2.16) 



in the sense of matrices and within Zpd{N, v) for sufficiently large iV, where 

j\{s)ds<Ci\og{e + t)y (2.17) 

JO 

for some constant C and all t > 0. 

If we assume (A3^), we also weaken (A4) to (A4^) 

Im Ait, + c|e|I + \l{t)\ + \5{t, 01 > (2.18) 

with the same 5(t,^) as in (A3^). 

Remark 8 This assumption is mainly for convenience and to simplify the 
main estimate in the pseudo-differential zone given in Lemma 9. If one wants 
to obtain sharp estimates in specific situations a more detailed analysis of the 
behaviour of solutions for small frequencies has to be carried out. We refer to 
[19] for an example, where such a more detailed consideration was essential. 

Another suitable replacement for (A4) could be to assume that the full symbol 
A(t,^) has all its eigenvalues in the (closed) complex upper half plane and that 
it is uniformly symmetrisable on Zpd{N' ,v). 

Later on we will discuss the construction of representations of solutions and re- 
sulting estimates under assumptions (Al) to (A4). Under the weaker assump- 
tions (A3") and (A4^) the uniform behaviour of the energy within Zj.eg{N, v) 
is replaced by an equivalence to exp Nq (5(r, Odrj and can be obtained by 
minor modifications in the arguments. 



3 Representation of solutions 



In order to solve (1) we apply a partial Fourier transform with respect to 
the spatial variables. Assuming Uq G iS(]R"), this gives a system of ordinary 
differential equations 

BtU = A{t,i)U, U{0,-) = Uo (3.1) 

parametrised by the frequency parameter ,^ G M" with initial data Uq G iS(M"). 
We construct its fundamental solution S{t, s, ^), i.e. the solution to the matrix- 
valued problem 

DtSit, s, = Ait, OSit, s, 0, Sis, s,0 = le C™^^'". (3.2) 

Then, the inverse partial Fourier transform gives f/(t, x) = Sit, 0, D)Uoix) for 
all t and x, while our construction gives a representation of £^(t, 0,D) as a 
matrix of Fourier integral operators. 



3.1 Treatment in the pseudo- differential zone 

We use the dissipativity assumption (A4) and restrict considerations to the 
zone Zpd{N,u). 

Lemma 9 Assume (A^)- Then the fundamental solution to (3.2) satisfies 

||£:(t,0,OII <Cexp(c"(log(e + t^))') (3.3) 

for some constants C and C and uniform in (t,^) G Zpd{N, v). 

PROOF. Fix ^ G M" and denote by V{t) the solution to T)tV{t) = A{t, OV'(t) 
with data V{0) = Vq. Then using the Euchdean inner product (■, ■) on C™ we 
obtain from (A4) that 

^^\\v{t)r = -2{imA{fov{t),v{t))<2cmv{t)r+imv{t)r (3.4) 

for all t < t^. Hence, by applying the Gronwall inequality the estimate 

\\V{t)f<\\Vofexp(2ct\^\ + j\{s)ds^<exp{C'{\og{e + t^)y) (3.5) 
follows. This proves the statement. 

Later on we need estimates for derivatives of S{t^,0,^) with respect to ^ as 
1^1 — )■ 0. They essentially follow from the symbol estimates satisfied by A{t,^) 
due to (Al) i^^i^. 

Lemma 10 Assume (Al)(_^^i^ and (A^)- For all multi-indices \a\ < min(£i,£2) 
the estimate 

l|Df^(t6 0,OII<C«ler'"'(log(e + t^))'^'"'exp(C"(log(e + t^))'^) (3.6) 

is valid uniformly on |^| < A^. 



PROOF. To prove this fact we use Duhamel's formula for ^-derivatives of 
(3.2). Let first |a| = 1. Then D^D^^ = (D^A)£ + A(D^^) and thus using 
D|'£(0,0,.^) = we obtain the representation 

D^£(t, 0, = f^Sit, T, 0(D^M(r, 0)^(r, 0, Odr. (3.7) 



Since ||D|'74(r, ,^)|| < 1, this implies the estimate 



< |er^(log(e + t5))%xp(c(log(e + t^ 



, 1/ 



uniformly on Zpd{N, v). 

For |q;| = £ > 1 we use Leibniz formula to represent D^£^(t, 0,^) by a cor- 
responding Duhamel integral using a sum of terms (D^M(t,,^))(D^^£^(t, 0,^)) 
for |q;i| + \a2\ < i, C(2 < a, and apply induction over i to obtain 

l|D5"^(t,0,e)||<|er'°'(log(e + t^))''"'exp(C"(log(e + t^))'), \^\ < N. 

(3.9) 

Derivatives of S with respect to t are estimated directly by the differential 
equation and ||A(t, ,^)|| < |,^|. Combined with 

iD^^d < k\^r < |er'-l"l (log(e + t^))' (3.10) 

this completes the proof. 



Remark 11 If one is only interested in energy estimates, the statement of 
Lemma 9 is sufficient to get the corresponding micro- energy estimate in the 
pseudo-differential zone. 

Remark 12 Any replacement of (A4), which is sufficient to deduce Lemma 9, 
is also enough to prove Lemma 10. This is interesting in particular for the 
treatment for second order scalar problems, where other ways to prove the 
estimate of Lemma 9 are available and important to obtain sharp result. See, 
e.g., the treatment in [19] or [4] where reformulations as integral equations 
were used. 



3.2 Treatment in the hyperbolic zone 



The following considerations are only applied within the hyperbolic zone. We 
can write formulae globally for all t and C, if we make use of the cut-off function 
Xhyp{t,C,). We omit this to keep notation as simple as possible. 



10 



3.2.1 Treatment of the main part 

We apply the diagonaliser M(t, of Ai(t, 0- For this we set U^°'> = M-\t, ^)U, 
so that 

Af/^°^ = {BtM-\t,0)U+M-\t,0^tU 

= [{DtM'')M + M-\t,OA{t,OM{t,0)u^''^ (3.11) 



is vahd. Now the main term V(t,C,) = diag(Ai(t,,^), . . . , \m(t,0) is diagonal 
and the remainder is given by 

Ro{t, = M'\t, 0{A{t, - A,{t, 0)M(t, + (D,M-i(t, 0)M(t, 0- (3.12) 



The rules of the symbohc calculus imply Ro{t,^) E S^^'J' {0,1}, while our 
assumptions on the characteristic roots give T){t,^) G S'^1^{1,0}. 

3.2.2 The diagonalisation scheme 

The initial diagonalisation step was done in the previous section. We trans- 
formed the system by a matrix M G (S^^^{0,0} to a diagonal form mod- 
ulo Sj^'i^'^^ {0)1}- III the k-th. step we want to diagonalise it modulo the 
class Sj^'i^'^^ ^ {—k,k + 1} for sufficiently large A^. (The choice of A^ de- 
pends on k in order to guarantee the invertibility of the transformation matrix 
Nk g5jv,40,0} below.) 

Recursively, we construct matrices Nk(t,^) to diagonalise the system Dt — 
V{t,^) — Ro{t,^) modulo iS^'j^^ {— ^j k + 1}. We use the notation 

N,{t,0 = ^N'^^\t,0, F,{t,0 = E^^'H^.O, (3.13) 

j=0 i=o 

where initially iV(°)(t,0 = I,5W(t,0 = i?o(^,0 and F(o)(t,0 = diag 5(o)(t,0. 
The construction goes along the following recursive scheme. 

Step k, k>l.We set 



F(^-i)(t,0 = diagi?(^-i)(t,0, (3.14) 

{ -B(fe-i)(f,g)ij . / . 

A,(t,0-A,(i,?)' '^^ (3.15) 

0, t=J, 

i?« (t, = (D, - V{t, - Ro{t, 0)Nk{t, 

- Nkit, 0(Di - Vit, - Ffc_i(t, 0). (3.16) 
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We prove by induction over k the estimates N^''\t,^) G S^^'^/ {—k,k} and 
B(-''\t, e S^^'f^'''''\-k, k + 1}. For A; = 1 we know 

i?«(t,0G5^V''~'{-i,2}, (3.17) 
the last one follows from the representation 

i?«(t,0 = D,iV«(t,0 - (i?o(t,0 - F^^\t,0)N^'Kt,0- 

For k > 1, assume we know B'^^^^^t,^) E Sj^'J' {—k + l^k}. Then by 
definition of iV(^) we have from (Ai(t,0 - Aj(t,0)"^ e '^^f {-1,0} that 
ArW(t,0 e 5^;'=^-'^{-A;, A;} and obviously F('=-^)(t,0 G 5^;'=^-'^{-A; + 1, A;}. 
Moreover, by algebraic manipulations we get 

k k k—\ 

5(^) = (D, - 1? - i?o)(E ^^'^) - (E ^'^'^)(Dt - ^ - E ^^'^^) 

k—l k 

+ A^^'') Yl ^^'^'^ - (E n'^''^)f'^''1 

u=0 u=l 

The definition of A^*^'^^ (t, ^) was chosen in such a way that we have 

B('-'\t,0 + [N^'\t,OMt,0] - F^'''\t,0 = 0. (3.18) 

The sum of the remaining terms belongs to the symbol class S^',/ {— A;, A; + 
1}. Hence we have proved B^''\t,^) G S^'j/ {—A;, k + 1}. Furthermore, the 
definition of B^''\t,^) implies the operator identity 

(d, - V{t, - Ro{t, 0)Nk{t, = Nk{t, (d, - V{t, - Ffe„i(t, O) (3.19) 

modulo 5^',^^" ~ {—k,k + 1}. Thus we have constructed the desired diago- 
naliser if we can show that the matrix N^it, ^) is invertible on Zhyp{N, v) with 
uniformly bounded inverse. But this follows from N^ — / G Sj^'J^~ { — 1,1} by 
the choice of a sufficiently large zone constant N . Indeed, we have 

l|iV.(t,O-/|l<C^^^^^^<^^0 asiV^oo. (3.20) 

Thus with the notation Rk{t,^) = —Ni7'^(t,^)B''''\t,^) we have proven the 
following lemma. 
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Lemma 13 Assume (Alji-^^i^ and (A2). For each 1 < k < £2 there exists a 
zone constant N and matrix valued symbols 

• Nk{t,^) G 5^'i^^~ {0,0}, which is invertible for all (t, ^) G Zhyp(N,iy) and 
its inverse satisfies N^^{t,^) G S^^'^/^ {0,0}; 

• Fk^i{t,^) G Sj^'J^~ {0, 1}, which is diagonal; 

. R,{t,OeSf^;',^~'-'{-k,k + l}; 

such that the (operator) identity 



Dt-vit,0-Roit,0)Nkit,o 

= Nk{t, (Bt - V{t, - Fk-i{t, - Rk{t, 0) (3.21) 



holds true for all (t,^) G Zhyp{N, u). 



3.2.3 Remarks on perfect diagonalisation 

Under assumption (Al) 00,00 Lemma 13 can be understood as a perfect diago- 
nalisation of the original system. If we define F{t, C,) as an asymptotic sum of 

theF('=)(t,0, 

00 

F{t,0-T.F^'\t^O, (3.22) 

k=0 

this means we require F(t,C,) — Fk(t,C^) G S]\[^^{—k — 1, /c + 2} for all k E N, 

and similarly 

00 

iV(t,0~E^^'ni,0, (3.23) 

A;=0 

which can be chosen to be invertible, equation (3.19) implies 

(d, - v{t, - R{t, 0)N{t, - N{t, (d, - v{t, - F{t, 0) e nN,M- 

(3.24) 
In this sense the symbol class 'Hn,u{^} is understood as the residual symbol 
class within the scheme of diagonalisation. 



3.2.4 Solving the diagonalised system in the regular sub-zone 

Our strategy is as follows. In the first step we solve the diagonal system D^ — 
V — Fk-i- Its fundamental solution £kit, s, ^) is given by 

Skit, s, = exp (i 1^ (Vir, + Ffc_i(r, 0)dr) 

= exp (^i^*I?(r,Odr) exp (^i^V,_i(r,Odr) . (3.25) 
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The second factor is uniformly bounded with uniformly bounded inverse on 
Zreg{N, z/) , bccausc Ffc_i - F(o) G 5^;,{-l, 2} C LfL]{Zreg{N, u)) by Propo- 
sition 2, (5), and ImF*^°) satisfies (A3). The first factor is a unitary matrix. 
Thus, we obtain the uniform bound 

^(t,S,0|| < C, (t,0,(s,0 e Zreg{N,u) (3.26) 

regardless of the order of s and t. 

In the second step we make the ansatz 

£k{t,s,0=£k{t,s,OQk{t,s,0 (3.27) 

for the fundamental solution Sk{t, s, ^) to Dt — T> — Fk — Rk- Then the matrix 
Qk(t, s,^) satisfies 

BtQkit, s, = T^kit, s, Qk{t, s, 0, Qk{s, s,0 = I^ C'"''™ (3.28) 

with coefficient matrix 

TZkit, s, = Skis, t, OMt, i%{t, s, 0. (3.29) 

Rk G SN,u{—k,k + 1} with k > 1 implies uniform integrability of TZk over 
the regular sub-zone. Therefore, the representation of Qk{t, s, C,) as the Peano- 
Baker series 

Qfc(t, s, = I + f; i' r 7^fc(tl, s, r 7^fe(t2, s, 
fc=l -^"^ -^^ 

•••y"'7^fe(4,s,0d4■■■dtl (3.30) 



implies boundedness of Qk(t,s,^) uniformly in {t,C,),{s,^) G Zreg{N,u). We 
obtain even a little bit more: 

Lemma 14 Assume (A1)q^^,^, (A2) and (A3). Then the matrices Qk{t,s,C,), 
k < £2, defined as solutions to (3.28) are uniformly bounded on Zhyp{N,v), 
invertible and tend locally uniformly in ^ ^ to Qk{oo, s, ^) G L°°{Zhyp{N, u)) 
as t —)■ 00. 



PROOF. The convergence follows directly by applying the Cauchy criterion 
to the representation (3.30). Invertibility is a consequence of the Liouville 
theorem, which gives 



det Qfe(t, s,^) = exp ( i / trace 7^fc(^, s,,^)d^ j = exp ( i / trace Rk{0,^)d9 
> exp (~m f \\Rki9,^)\\d9] > C > 0, (3.31 
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uniform in (t,^), (s,^) G Zhyp{N, u). 

The following theorem explains the true nature of assumption (A3) provided 
that z/ = 0. 

Theorem 15 Assume (AI)o^od with z/ = and (A2). Then (A3) is equivalent 
to the existence of a constant C > such that 

C'^\\V\\ < \\S{t,s,^)V\\ < C\\V\\ (3.32) 

for all vectors V e C™ and (t,0, {s,0 e Zhyp{N,0). 



PROOF. Lemma 14 in combination with (3.26) gives the uniform bound 
under (A3). Without (A3) the estimate of (3.26) has to be replaced by a 
polynomial bound 

\\g,{t,s,Ol\\^k{s,t,0\\<Ck[^^ , t>s, (3.33) 

where the constant K is independent of k (and depends only on the first 
diagonal lower order term F°). Similarly, we obtain with the same K 



1+t 



K 



\\S,it,s,m < exp ^* II Im(Ffc_i(r,0 + i?fc(r,0)||dr) < C, Q-^ ^ 

'(3.34) 

for all t > s. Choosing k large, the polynomial decay of the remainder Rk{t, ^) 
becomes strong enough to compensate the increasing terms and we obtain the 
Duhamel representation 

Skit, s, = Skit, s, ozki-s, O-^l Skit, e, ORkie, O^kie, s, Ode (3.35) 
Zkis, = I + i / Skis, e, ORkie, O^kie, s, O^e < 1. (3.36) 

Indeed, the integral in (3.35) is convergent for k > 2K, 

) 

Skit,e,ORkio,OSki0,s,Ode 






and bounded by |^|"'^(1 + t)^-'=(l + s)"^. Similarly \\Zkis,Q-l\\ < |Cr'=(l + 
s)~^ and hence the first term has the lower norm bound ||£^fc(t, s,^)|| > (1 + 
t)^(l + s)~^ for fixed |^|. Choosing s big enough depending on |^| implies 
that Skit, s, C,) is a small perturbation of Skit, s, C,). 
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Assume now that (A3) is violated. Then we find sequences t^ — > oo, s^, and ^^ 
such that one matrix entry of the integral in (2.12) tends to either oo or — oo. 
We consider the +00 case, and assume w.l.o.g. that s^ > s for sufficiently 
big s and that the matrix entry corresponds to the first diagonal element. 
Then with ei the first basis vector £^fc(t^, s^,^^)ei — ?> 00 and therefore also 
S(t^, s^,^^)ei ~ Sk(t^, s^,^^)Nk{s^,^^)M{s^,^^)ei -> 00 which contradicts to 
the uniform upper bound. Similarly, the —00 case contradicts to the lower 
bound and the statement is proven. 

Remark 16 A similar argument does not apply if u > 0. In this situation the 
polynomial bound becomes a superpolynomial bound by exp(C(logt)^"^'^) which 
cannot be compensated within the diagonalisation hierarchy. 

The following results explain why we performed k > 1 steps of diagonalisation. 
In order to prove dispersive estimates later on, we have to control derivatives of 
Qk(t,x,C,). Formally differentiating (3.30) includes ^-derivatives oi Sk(t,s,^), 
which are increasing in t. This can be compensated by better estimates on 
the remainder Rk{t, C,)- As we will see below, more steps of the diagonalisation 
hierarchy allow for symbol-like estimates for more derivatives of the entries of 

Lemma 17 Assume (Al)^^^^^, (A2) and (A3). Then the matrix Qk{t,s,^) 
with 2k < £2, satisfies for all \a\ < min(A; — l,£i) 

D^"Q.(t,s,0||<ler'"'(log(e+|er^))'°' (3.37) 

uniformly in s,t > i^. Furthermore, for \^\ < N and \a\ < min(^^,£i), we 
have 

uniformly in t > t^. 



D?Q.(t,t5,0 <lerl°l(log(e+|er^))'"' (3.38) 



PROOF. The statement follows similar to the considerations in [10, Lemma 
2.10]. We concentrate on the second estimate, the first one is analogous. Note 
first, that by differentiating (3.25) we obtain 

mSk{t,k,0\\<{l + tpUg{e + tW"[ (3.39) 



the logarithmic term arises from differentiating Fk-i{t,C,) in (3.25). This im- 
plies that TZk{t, te.j satisfies weaker estimates than Rk(t, C,), for which we have 
Sj^' ^ {—k,k + 1} estimates. Based on the equivalent to [19, Proposition 
11] we obtain 



a\ 



|D^"7^,(t,t^,OI| < |er'-'"l(l + t)-^(log(e + t))'"' (3.40) 
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for \a\ < ^Y' ^"^^ W\ — ^1- Differentiating the series representation (3.30) 
term by term, taking into account these estimates and combining them with 

|D^"td^^~der"<ler'-l"l(log(e + t^))''^ (3.41) 

proves the desired resuh. 

Remark 18 If we diagonalise perfectly modulo 'HN,iy{i} we can construct the 
corresponding fundamental solution Soo(t, s,^) = £oo{t, s, C,)Qoo(t, s, ^) and the 
previous lemma holds for all multi-indices a. Note, that T-Ln,u{^} is invariant 
under multiplications by Sodt, s, ^). 



3.2.5 Estimate in the oscillating subzone 

It remains to estimate the fundamental solution in the oscillating sub- zone. In 
this part of the phase space it suffices to apply one step of diagonalisation, i.e. 
diagonalisation modulo 5^',^^" { — 1,2}. We construct the fundamental solu- 
tion £i{t,s,C,) = Si(t,s,C,)Qi (t, s,^) oi Dt — V — Fi — Ri following the lines of 
the previous subsection replacing the uniform integrability of the remainder 
by the estimate 

rk .._ rk (log(e + t)j 

iTtm 

<N[\og{e + t^)y. (3.42) 

k 



(log(e + t) 



This implies 



\\Qiit,s,0\\<ew(C\\ogie + t^)) ) (3.43) 

for t^ < s <t <t^. For later use we need estimates for derivatives of Qi with 
respect to ^. 

Lemma 19 Assume (Alj^^^^f^, (A2) and (A3), with £2 > 2. Then the matrix 
Qi{t, s,^) satisfies 



\\B^^Qi{t,t^,0\\ < \^r^''K^og{e + t)f^"'^^^^exp[C'[log{e + t^)) ) (3.44) 
for all \^\ < N and all multi-indices \a\ < ii. 



PROOF. See, e.g., [10, Lemma 2.9]. The proof follows directly by differenti- 
ating (3.30) (for k = 1), using the symbol estimates together with (3.10) and 
estimates (3.42). 
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3. 3 Combination of results 



We combine the representations of the two previous sections. We distinguish 
between three cases. If |.^| > A^ then only the hyperbohc zone occurs. Hence 

S{t, 0, = M-\t, ONk\t, 0^fc(t, 0, OiVfc(0, OM(0, 0- (3.45a) 

For |.^| < A^ and t <t^ only the pseudo-differential zone is of interest and the 
fundamental solution is constructed in Section 3.1. For t > t^ we obtain 

S{t, 0, = M-\t, ONk\t, OSk{t, t^, ONkit^, OM{t^, O^m, 0, (3.45b) 

with k = 1 in the oscillating sub-zone and large k and t^ in place of t^ in the 
regular sub- zone. 



We will bring these representations into a unified form. For this, let us define 

7' 

/o 
Then the following statement holds true. 



^.(t,0 = y^*A,(r,Odr. (3.46) 



Theorem 20 Assume (Al)i^^i^, (A2), (A3) and (A^)- Then the fundamental 
matrix S{t,0,^) can be represented as 

m 

^(t,0,O = Ee"'^'^^'^*'^^'^'^^.(^'O, (3.47) 

with matrices Bj(t,^) G C"^™ subject to the estimates 

\\B,{t,0\\ < Cexp(C"(log(e + t))') (3.48) 

tn Zpd{N, v) U Zosc{N, v), and 

D^fi,(t,0||<a|er'"'(log(e + t))(2^+^)Hexp(c"(log(e + t))'^) (3.49) 
for all \a\ < min(£i,£2/2 — 1) in Zreg{N, v). 



PROOF. Note first, that for (t,0 e Zpd{N,u) U Zosc{N,u), all the terms 
g±it|CI0j(t:?/l?l) are uniformly bounded and satisfy the corresponding symbolic 
estimates. Thus artificially introducing these terms does not destroy our state- 
ment. 

The estimate in the pseudo-differential zone and the oscillating sub- zone fol- 
lows directly from Lemma 9 and equations (3.26) and (3.43). Note that for 
z/ = we obtain just uniform bounds. 



In the regular sub- zone the matrices Bj collect terms from the diagonalisers 
M, iVfc, the matrices S{t^,0,C,) and Qk{t,C,) and the terms arising from 

expfif Ffc(r,Odr), (3.50) 

y Jmax{tj,0) / 

all of which satisfy symbolic estimates due to Lemma 10, 17 and for the 
last one due to F^ E S^^'^^^' {0, 1}. Note further, that (3.41) implies from 

||D^"7Vfc(t^,e)|| <C|erl"l(log(e + t5))''^'"' (3.51) 

which gives the 2v log-estimate in the regular zone despite the better estimate 
of Lemma 3.41. 



4 Examples 



We will give some examples to illustrate the applicability of our previous 
construction. 



4-1 Differential hyperbolic systems 



If we consider differential systems, they can always be written as D^t/ = 
A(t,D)f/, with 

n 

A{t,D) = Y.A,{t)J}^^+B{t), (4.1) 

i=i 
for t-dependent matrices Aj{t), B(t) G C™'^'". Assumption (Al)^^^^^ is satisfied 
if Aj{t) e %{0} and B{t) e %{!}, where 

TAp} = |/ e C-(M+) : \d^f{t)\ < C, ^-L-[log{e + t)yy '| , (4.2) 
together with the assumption that the associated homogeneous matrix 

n 

E^.WO (4.3) 

i=i 

has real eigenvalues for all t and ^. The system is symmetric hyperbolic if 
the latter matrix is self-adjoint, which means that all matrices Aj{t) are self- 
adjoint. We denote the eigenvalues of (4.3) by Ai(t, ^), . . . , \m(t, ^) in ascend- 
ing order. Assumption (A2) guarantees that they are uniformly distinct if 
restricted to M x S". 
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For the following we assume the system to be symmetric hyperbolic. Then 
Assumption (A3) simplifies to 



sup 

s.t 



diag M-\e, 0(lm B{e))M{e, d^ 



< CXD 



(4.4) 



for the (0-homogeneous) unitary diagonaliser M{t, ^) of (4.3), while (A4) reads 
as 

Im5(t) +7(t)/> (4.5) 

for a suitable function 7(t) satisfying (2.17). Both assumptions are trivially 
satisfied if we assume lm.B{t) G L-'^(]R+). 

This corresponds to the situation in [3] . Main difference between the results of 
[3] and our results is, that we do not subsequently require the non-vanishing 
of the Gaussian curvature of the characteristics. In Theorems 26 and 27 we 
derive estimates for solutions without this extra assumption. The case of [3] 
is covered by the case 7 = 2 in Theorem 26. 



4-2 Hyperbolic equations of second order 



Second order equations under related assumptions are considered in [10], [11] 
and related papers. If we consider the equation 

n n 

B^.u = DtJ2D,^^bj{t)u+ J2 aij{t)B^fi.^M+c{t)BtU+J2djit)I)^M + e{t)u 
j=i i<*<i<" .7=1 

(4.6) 
with aij(t),bj{t) G 71,(0}, c(t),dj(t) G %{!} and e(t) G Tu{2}, we can rewrite 
it as 2 X 2 pseudo-differential system in f/ = {h(t, D)-u, D^-u)-^, h(t, ^) ~ |.^|i (be- 
ing ^-independent within the pseudo-differential zone and smoothly continued 
to the hyperbolic zone by adding a suitable term from 'Hn,u{^}), 



DtU = A{t,D)U, 



(4.7) 



where within Zhyp{N, v) we have 



A{t,0 



lei 

je|-'Eo.,,(t)e4E^.We,, 





, 

+ 1 (4.^ 

\i\-Ht) 0. 
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modulo lower order terms form 'HN,iy{^}- The second line gives the terms from 

<Sn,u{0, 1} and Siy^,,{—1, 2}, respectively. Similarly, we obtain within Zpd{,N, v) 

hit) 

I DtXoghit) \ f 0\ 

+ + (4-9) 

\h{t)-^Y.d,m,c{t)] \h{t)-^e{t)Q^ 

with h{t) = h(t,0) = YT^(log(e + t)y. Note that some of the 'lower order 
entries' are the dominant terms now. 

Assumptions (Al)£^^^2 ^'^d (A2) are satisfied if we require that eigenvalues of 
the homogeneous principal part 



1 " 



\ 



Ujlbjmj) +E«MWe.e, (4.10) 



are real and distinct, i.e., if bj{t) is real and 

0<(j2bjmj) +4E«MW6e,, ^7^0. (4.11) 

Assumption (A3) requires that the integral of 

^^'^^ imc(t)+ ^^ '^; + r^' ; ') (4.12) 



A+(t,0-A-(i,0 "' A+(t,0-A-(t,0 A+(t,0-A_(t,0 

is uniformly bounded over the hyperbolic zone. This is in particular the case 
if Imc(t) and Im(ij(t) are integrable and the equivalent to assumption (4.16) 
is satisfied. Furthermore, (A4) becomes a condition on the lower order terms 
guaranteeing that their influence is dominated by the principle part. 



J^.3 Hyperbolic equations of higher order 

As third example we consider hyperbolic equations of higher order, 

m— 1 

D>+E E a„,(t)Dp> = 0, (4.13) 

i=0 j+\a\<m 

together with the Cauchy data 

Dlu{0,-)=u„ j=0,l,...,m-l. (4.14) 
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Similar to the second order case, this can be brought into pseudo-differential 
system form. 

Then {Al)^-^^^,^ is satisfied if we assume aj^^ £ Tlim — j — \a\} and if the roots 
r = \j{t,^) of the homogeneous principal part, i.e., the solutions of 

m— 1 

r™+E E aatyC = 0, (4.15) 

i=0 j-\-\a\=m 

are real. For (A2) we require that they are distinct. If we consider equations 
which are homogeneous of order m, i.e., aj^ait) = ioi j + \a\ < m, condition 
(A3) is equivalent to 



max sup sup 






< oo. (4.16) 



while Assumption (A4) is trivially satisfied. 



Remark 21 In [6] homogeneous equations of order m where considered un- 
der the stronger assumption dtaj^ait) G L"'^(]R+) for j + |a| = m. This implies, 
in particular, that (4.16) is satisfied. The representation of solutions is ob- 
tained there by the asymptotic integration method, requiring less conditions on 
the coefficients aj^a (conditions on only the first time derivatives are enough). 
However, the asymptotic integration method yields a representation of solu- 
tions with less control on its amplitudes and an additional loss of regularity in 
the dispersive estimates, which does not occur with the present method. 



5 Dispersive estimates 

In this section we are concerned with dispersive estimate for solutions repre- 
sented in the form of Theorem 20. From now on we assume that n > 2. 

First, we give the estimate for low frequencies: 

Lemma 22 Assume (A4). Then solution U = U{t,x) to (1) satisfies 

<C(l+t)-"(log(e + t))'"%xp(C'log(e + t)'^)||f/o||Li(R") (5.1) 
localised to Zpd{N, v) U Zosc{N., v) (for any choice of N). 

PROOF. Based on the mapping property of the Fourier transform J-" : L^(R") — )■ 



22 



L°°(M") and Holder inequality it is sufficient to estimate 

ll^(^) 0,OXpd+osc{t,q)\\L'^(M.r^) < \\S{t, 0,.^) II ^oo(|g|<|^) II Xpd+osc II Li(R") 

and, therefore, the estimate follows from Lemma 9 and the geometry of the 
zone. 



For the high-frequency part we at first derive some abstract statements giving 
LF-L'^ decay rates for oscillatory integrals with a related structure to the ones 
constructed in Theorem 20, 



Ttf{x)= I e'(^-«+*^(*'«6(t,0/(0de, 



with a real phase function ip and and amplitude h. For simplicity, we omit the 
inclusion of logarithmic terms in assumptions here, later they will give a (for 
V < 1 small) change in the decay rates by a simple argument. 

We introduce a cut-off function of the form '?/'((l -|- |t|)0 ^r some ip G C^{W^) 
such that ip{i) = 1 for |,^| < i, and for |,^| > 1. We recall that we use the 
notation L^(R") for the homogeneous Sobolev space 1V^(M"). 

For now we may assume that phase functions satisfy ip{t, ^ > uniformly for 
all ^ 7^ in the support of b{t, C,)- This can always be achieved by localisation 
and in combination with adding a linear function to the phase which will not 
affect estimates. Dispersive estimates for the corresponding Fourier integral 
operator Tf are determined by the geometric properties of the family of Fresnel 
surfaces 

S* = { e G M" \ {0} : ^(t, = \^Mt, e/lel) = 1 }. (5.2) 

Following [18] we introduce two indices for such a Fresnel surface S, assuming 
that it is of class C^ with k being sufficiently large. For -u G S we denote by 
Tu the tangent hyperplane to S at u. Then for any plane H of dimension 2 
which contains u and the normal of S at m we denote by 7(S; u, H) the order 
of contact between the curve S fl // and its tangent H ClTu- Furthermore, we 
set 

7(S) = supsup7(S; -u, iJ), 7o(S) = supinf 7(S; -u, iJ). (5.3) 

u H u H 

Obviously the definition implies 2 < 7o(S) < 7(S). For isotropic problems S 
is a dilation of E>"'~^ and 'j{E>"'~^) = 7o(S"~"'^) = 2. Moreover, if the Gaussian 
curvature of S never vanishes we have 7(S) = 2. 

To control the order of contact of E by tangent lines quantitatively, we will 
introduce another quantity x(S). First, assume that S is convex and that it is 
of class C"'^^'^'^^. For u G T,, rotating E if necessary, we may assume that it is 
parameterised by points {{y, h{y)),y G fi} near u for an open set fl C M"^^. 
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For u = {y, h{y)), let us define 



7(S) 

x(S; u) = inf ^ 



1^1=1 



i=2 



— /i(y + pcj)|p=o 



From the definition of 7(S) it follows that >f(S; -u) > for all u ^ Q. Indeed, 
from the definition of 7(S; u, H) it follows that if u is such that y + pu E H, 
then 

^7(S;n,//) 



x(S; -u, iJ) 



>0. 



Now, we clearly have I]J=2 -§-7h{y + p(^)\p=o > x{Tj;u, H), and hence we 
have >i{T,;u) > since the set |a;| = 1 is compact. Noticing that x[Tj;u) is a 
continuous function of m, by compactness of S it follows that if we define 



x(S) = min x(E; -u) 



then xfE) > 0. 



If S is not convex, we define 



xo(S) 



70 (S) 

min sup V" 



dpi 



h{y + pu)\p=o 



Again, we have xo(S) > 0. 

The quantities x(S*) and xo(S*) evaluated for all times t will allow us to ensure 
that if the surfaces degenerate, it will be at least uniformly with respect to 
t. If the degeneracy of a certain order is not uniform, it may become uniform 
when we increase integers 7 in the formulations below. 

For completeness we include the following proposition. It gives a criterion on 
the convexity of the level sets of the phases, as well as an upper bound on the 
index 7 for systems arising from differential operators. 

Proposition 23 Assume that detAi{t,C,) is polynomial in C, and denote by 
Xk{t,^), k = 1, ... ,171, the characteristic roots of the operator (1), ordered by 
Ai(t, ^) < A2(t, .^) < ■■■ < \m(t,^) for ^ 7^ 0. Suppose that all the Hessians 
V|Afc(t, ^) are semi-definite for ^ 7^ 0. Then there exists a polynomial a(t, {) 
in ^, of order one, such that \m/2(f,0 < ^(t,^ < ^m/2+i(^)0 {if 'm is even) 
or a{t,^) = A(m+i)/2(^)0 {if ITT- is odd). Moreover^ , the hypersurfaces S^ = 
(e e M"; Mt, = ±1} yj^th ^k{t, = Mt, - a{t, {k^{m + l)/2) are 
convex and we have 7(S|,) < 2[m/2j. 



Here, as in (3.46), we define (pj{t,^) = j Jq Aj(t, ^) dr. 
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By taking the determinant of the system (1) we can reduce the analysis of 
characteristics to scalar equations. Then this proposition can be readily shown 
by a modification of the argument in [17], so we omit the details. 

We are now ready to estimate the appearing oscillatory integrals. In the sequel, 
for r > 0, by \r\ we denote its integer part. The Proposition below extends 
Proposition 4.2 in [6] where the limit of (f{t, ^) was assumed to exist as t — ?■ oo, 
and which deals with amplitudes of type (0, 0). 

Proposition 24 Let 7 G N. Let Tt be an operator defined by 

TJ{x) = f e^(-«+*^(*'«)) [1 - ^((1 + |t|)0] b{t, m de, 

where ip{t, C,) is real valued, continuous in t, smooth in ^ E M"\0, homogeneous 
of order one in C,, and such that for some to > and C > we have 

C-^lei < V{t,0 < C\^\ and |9^Xt,0l < C|er-I°l for all t > to, C 7^ 0, 

(5.4) 
and all a such that \a\ < max{7 + 1, [(n — l)/7j + 2}. Assume that the sets 

S* = {eGM"\0:v^(t,O = l} (5.5) 

are convex for all t > to, and assume that hmsupj^oo7(S*) < 7 and that 
liminfj_^oo ^(S*) > 0. Let us suppose that the amplitude b{t,^) satisfies 

\dpit,0\<C^\^r^''^ for all |a| < [(n - l)/7j + 1. (5.6) 

Let l<p<2<q<oo be such that ^ + ^ = 1 . Then for all t > to we have 
the estimate 

n-l /I 1\ 
\\Ttf\\Ll{R") < Ct ^ yP ^^||/||i^^(Kn), 

where Np = (n - ^) (^ - j) • 



PROOF. Since the L^-estimate ||Tt/||i2 < C||/||l2 readily follows from the 
Plancherel identity, by interpolation we only need to prove that 

— "~^ 

||^t/||L°°(R") < Ct -y ||/||li(R"), (5.7) 

with Tt = Tf o |Z)|^^i, where the amplitude a{t,^) of Tt satisfies \d?a{t,C,)\ < 
C«|^|-^i-l"l for all |a| < [(n - l)/7j + 1 and A^i = n - ^. Let us make an 
additional decomposition 

ft = Tt^'^ + t/') = ft o{1- ^{D)) + Tt o ^{D). (5.8) 

First we will treat the operator T/ . To simplify the notation, for this part we 
will denote a(t, ^)(1 — ip{^)) by a{t,^) again. By using Besov spaces, we can 
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microlocalise the desired estimate (5.7) to spherical layers in the frequency 
space. Indeed, let {$j}°2=o ^^ ^^^ Littlewood-Paley partition of unity, and let 

(oo 

be the norm of the Besov space B^^. Then, because of the continuous em- 
beddings LP C 5° 2 for 1 < p <' 2, and 5° 2 C L" for 2 < g < +00 
(see [1]), it is sufficient to prove the uniform estimate for the operators with 
amplitudes a{t,^)^j{^). Writing $j(0 = ^jiO^i^^^) with some function 
\E' G C^(0, 00), we may prove the uniform estimate for operators with ampli- 
tudes a(t, ^)^( 23 )• Such choice of \1' is possible due to our assumption (5.4) 
on ip(t,^), and we restrict the analysis for large enough t. Let 

I{t,x)= f e^(-«+*-(*'«))a(t,0^f^V^ 

JR" \ I-' J 

be the kernel of the corresponding operator. Since we easily have the L'^-L^ 
estimate by the Plancherel identity, by analytic interpolation we only need 
to prove the L^-L°° case of (5.7). In turn, this follows from the estimate 

_ n— 1 

|/(t,a;)| < Ct -I , with constant C independent of j. 

Let K, G C^(]R"') be supported in a ball centred at the origin, with some radius 
r > to be chosen later. We decompose the kernel I{t,x) as 

I{t,x) = Ii{t,x) + l2{t,x) (5.9) 

= /^„ e^(-«+*-(*'«)a(t, O'^ {t-'x + VMt, 0) ^ (^) de 

+ l^ e^(-^+*-(*'^))a(t, (1 - «:) {f-'x + V,v^(t, O) ^ (^) d^ 

n-l 

We can easily see that |/2(t,x)| < Ct t . In fact, we can show |J2(t, a;)| < 
Ct~' for / = \_{n — l)/'y\ +1 and then the required estimate simply follows since 
/> (n—l)/7. Indeed, on the support of l — K, we have |x+tV^V9(t,.^) I > rt > 0. 
Thus, integrating by parts with operator P = ^i^i^v^ (f f)p ' ^?' ^^ S^t 

h{t,x) = 

V(i,0^ 



J^y(.<+Mt,o)^P*y ^(^^^)(i _ ^) (t'^x + VMt,0)'^\ 



2^ 



d^. 
(5.10) 



Using the fact that \d^(f{t,^)\ < C\^\^ '"' by (5.4), we readily observe from 
(5.10) that the required estimate |/2(t,x)| < Ct~^ holds. Here we also used 
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the condition (5.6) assuring that we can perform the integration by parts 
L(n — l)/7j + 1 times. We note that since there is one more V^ip involved 
here, the condition \a\ < [(n — l)/7j + 2 in (5.4) allows us to integrate by 
parts / times. 

Now we will turn to estimating /i(t, x). Here we are going to use the structure 
of the sets S* in (5.5), restricting to t large enough. We recall that (5.4) 
implies, in particular, that f{t,^) > for all ^ 7^ 0. By rotation, we can 
always microlocalise in some narrow cone around e„ = (0,...,0,1), and in 
this cone we can parameterise 



{{y,ht{y)):yeU} 



for some open U C M"^^. In other words, we have ip{t; y , ht{y)) = 1, and it 
follows that hf is smooth and Vht : U — ?■ Vht{U) C W^~^ is a homeomorphism. 
The function ht is concave if S* is convex. We claim that 



I'^yhtiv)] < Ca, for all y E U and large enough t, (5.11) 



for all \a\ < max{7 + 1, [(n — l)/7j + 2}. Indeed, let us look at \a\ = 1 first. 
From ip(t; y, ht{y)) = 1 we get Vyip + d^^ip ■ Vht{y) = 0. From (5.4) we have 
|V^9?| < C, so also \Vyip\ < C. By Euler's identity we have d^^ip(t;en) = 
ip{t\ e„) > 0, so we have l^^nV^I > c > since we are in a narrow cone around 
e„. From this it follows that \Vyht{ii)\ < C for a\\y eU and t large enough. A 
similar argument proves the boundedness of higher order derivatives in (5.11). 

Now, let us turn to analyse the structure of the sets S*. We have the Gauss 
map 

|VcV'(t;C)| 

and for x = {x', a;„) G W''^'^ x M near the point — V(^v9(t; e„) we define Zf E U hj 
{zt,ht{zt)) = G''^{-x/\x\). Then {-Vyht{y),l) is normal to S* at {y,ht{y)), 
so we get 



_x_ _ {-Vyhtjzt),!) x^ 

\X\~ \i-Vyhtiz,),l)\ '"'' Xr.~ ''' 



and — = -VyhAzt). 



Making change of variables ^ = {\y, \ht{y)) and using {p(t, ^) = A, we get 
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h(t,x) 



e^M^'-y+-^My)+t)a^t^Xy,\htiy))^{^]K,it,x,y) 



JU 

oo 



JU 



JX(~XnVyht(zt)-y+Xnht{y)+t) 



X 



2J 



\^a{t,\y,\ht{y)) 



d^ 



d{X,y) 



dy dX 



X 



X ^ ( ^ J /«o(t, a;, y)xit, y) dy dX 



JU 



x^l 



^^X(-^yh^iz,)■y+h,{y)+tx-')^^^ ^^^ ^y ^ A/li(y)) A""^-' X 



2^X„ 



\x. 



-n+l+l-l 



fio{t,x,y)x{t,y)dy dX, 



where K,o{t, x,y) = n (t ^x + V^ip{t; y, ht{y))), and 

a{t,Xn,Xy,Xht{y)) = {x'^XYa (t,x;^^Xy,x;^^Xht{y) 



and where we made a change A = x„ A in the last equahty. Here also we used 



d? 



= A" ^xitjV)} where xi^iU) ^"^^ its derivatives with respect to y are 
bounded because of (5.11). 

If we choose r in the definition of the cut-off function k sufficiently small, then 
on its support we have \x\ ~ \xn\ ~ t, and we can estimate 



\h{t,x)\ <Ct 



-n+l 



^(A...)*(^)a— ' 



dA 



Qj^-n+l2Jin-l) 



X 



J{2^X,z,)<^\-\X 



n-l-l 



(5.12) 



dA, 



with 

J{X,zt) = |^e'^(-^-'^*(^*)-^+'^*(^)+*^")a(t,a;„. Ay, A/i<(t/))«:o(t,a;,i/)x(t,i/) dy. 
We will show that 

n-l 

\J{X,zt)\<C{l + Xy—, A>0. (5.13) 

Then, if we take I = n — ^^— i, and use (5.12) and (5.13), we get 



n — 1 ■ Ti — 1 



n— 1 / A \ n—1 1 

-< I (2-'A)"— ^ - A— ^MA 



\h{t,x)\ <cr—2^- 

= Ct~— X-^^i-\dX = Ct-— A-^^(A)dA 



(5.14) 



<cr 
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which is the desired estimate for Ii{t,x). 

Let us now prove (5.13). We note first that with this choice of / we have 

\dya\ < C for all |a| < [(n - l)h\ + 1. (5.15) 

Now, estimate (5.13) follows from Theorem 29 in Appendix A. Indeed, we can 
write J(A, Zt) in polar coordinates (p, cu) with y = pu + Zt, so that 

J(A, 2i) = e'^(*"" +'^*(^'» / / e'^^'^'''''''^f3(p,zt,u)p''~^ dpdu, (5.16) 

J§"-2 JO 

with 

F(p, Zt, u) = ht{pu + Zt) - ht{zt) - pVyht{zt) ■ u, (5.17) 

/3{p, Zt, u) =a {t, Xn, \{puj + Zt), \ht{pu + zt)) K,o(t, X, pu + zt)x(t, pu + Zt), 

where we can assume in addition that x = unless pu + Zt G U, so both p 
and u vary over bounded sets. Now, we want to apply Theorem 29 to obtain 
estimate (5.13). The function F in (5.17) satisfies condition (F3) of Theorem 
29 because of the definition of the convex index 7 and because ht is concave. 
From the assumption liminfj^oo ^iX't) > it follows that function F satisfies 
property (F2) of Theorem 29. One can readily see that the other conditions 
of Theorem 29 are satisfied, implying (5.13). 

We now prove the estimates for the integral T^ in the decomposition (5.8). 

This proof is essentially similar to that for T/ with a few differences that we 
will point out here. Again, by interpolation, it is sufficient to prove estimate 

llTi'VllLoo(Rn) < Cr'T^i|/|Ul(Mn), (5.18) 

with amplitude a{t, ^) satisfying 

|a^"a(t,OI<C,|er^^~'"' for all |a| < [(n - l)/7j + 1, (5.19) 

with Ni = n-^^^. 

We continue to estimate the integrals Ii{t, x) and l2{t, x) in the decomposition 

(5.9) corresponding to the integral T^ . In general, since for Tj we work with 
low frequencies |^| < 1 only, no Besov space decomposition is necessary, so we 
do not need to introduce function \E' and $j, and we can take \E' = 1. 

The additional complications are related to the fact that in principle deriva- 

fives of the amplitude of the operator T/ may introduce an additional growth 
with respect to t. In the estimate for /2(t, x) we performed integration by parts 
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with operator P. Now after the integration by parts the amphtude of this in- 
tegral in (5.10) is 



(P 



*\i 



[1 - v^((i + itmw, 0(1 - k) {t-'x + vMt, 0) ^ ( ^ 



Then, if any of the ^-derivatives falls on [1 — ip{{l + 1^1)0]) "^^ S^t an extra 
factor t which is cancelled with t~^ in the definition of P. However, in this case 
we can then restrict to the support of Vip which is contained in the ball with 
radius (1 + |t|)^^, so we are in the situation of very low frequencies \C,\ < Ct^^ 
again. Consequently, for this integral we actually get a better decay rate of 
Lemma 22. If none of the derivatives in (P*)' fall on [1 — ■?/'((l + jtDO], the 
argument is the same as in the proof of the estimate for l2{t, x) for the integral 

(2) 

Let us now analyse the term Ii{t,x) corresponding to T^ . Recall now that 

in the process of estimating /i(t, x) corresponding to T^ , we made a change 
of variables A = x^^X. As it was then pointed out, if r in the definition 
of the cut-off function k is chosen sufficiently small, on its support we have 
l^nl ~ |t|. On the other hand, we have \C,\ ~ A by the definition of A, since 
we may assume that fit,^) is strictly positive for ^ 7^ 0. It then follows that 
(1 + |t|).^ ~ A|x„| ~ A, and so the change of variables A = x~^\ changes 
[1 — ip{{l + 1^1)0] i'^to 1 — ip{\) in the amplitude of Ii(t,x). Justifying this 

argument, we can then continue as in the case of T^ . The crucial condition for 
the use of Theorem 29 is the boundedness of derivatives of a in (5.15). Here, 
every differentiation of a with respect to y introduces a factor x~^X which is 
then cancelled in view of assumption (5.19). It follows that \_{n — l)/7j + 1 
y-derivatives of a are bounded, implying the conclusion of Theorem 29. This 
yields estimate (5.18) in the way that is similar to the proof of the same 

estimate for T} . 



We now give an analogue of Proposition 24 without the convexity assumption. 
Naturally, in this case we get slower decay rates at infinity. 

Proposition 25 Let 70 G N. Let Tf be an operator defined by 



where ip{t, ^) is real valued, continuous in t, smooth in ^ E M"\0, homogeneous 
of order one in ^, and such that for some t^ > and C > we have 

C~'\^\ < ^(t,0 < C\^\ and |9fv^(t,0l < C|er-I"l for all t > to, ^ 7^ 0, 

(5.20) 
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and all a such that \a\ < 70 + 1- Assume that the sets 

S* = {eGM"\0:v^(t,O = l} (5.21) 

satisfy limsup^__^oo7o(S*) < 70 and lim inf t_^oo >^o(S*) > 0. Let us suppose that 
the amplitude b{t, C,) satisfies 

\dp{t,0\ < Caier'"' for all \a\ < 1. (5.22) 

Let l<p<2<q<oobe such that ^ + ^ = 1. Then for all t > t^ we have 
the estimate 



^hereN,= {n-^)[l-^^. 



PROOF. Let us show how the proof of Proposition 25 differs from the proof 
of Proposition 24. First, we need to prove that \I{t,x)\ < Ct ''o , t > to, for 
Ii(t,x) as in (5.9). We note that 70 + 1 > 2 (actually, as we observed before, 
we must have 70 > 2), so to prove the estimate for l2{t,x) we can show that 
\h{t,x)\ < Ct^^. This can be done by integrating by parts with the same 
operator P one time, and using (5.22) instead of (5.6). As for the proof of the 
estimate for Ii{t,x), we can reason in the same way as in Proposition 24 to 
arrive at the estimate (5.12), i.e.. 



\h{t,x)\ < c'r"+'2-''("-') 

/o 



J(2^A,^,)vp(^jA"-i-' 



d\ 



with the same operator 



J{\,zt) = / e'^^-^y'^^^^^>y+^^^^y^+'^-\{t,Xn,\yAhtiy))no{t,x,y)x{t,y)dy. 
Ju 

Now, instead of (5.13) we will show that 

\J{X,zt)\<C{l + X)~^, A>0. (5.23) 

i_ 

Then, taking / = n — -, we get the estimate |/i(t, a;)| < Ct 10 in the same way 

as in estimate (5.14). Now, estimate (5.23) follows from Theorem 29 in the 

appendix with N = 1. Indeed, let us write J{\,Zt) in the form (5.16)-(5.17) 

with phase 

F{p, zt, oj) = ht{pu + zt) - ht{zt) - pVzht{zt) ■ u. 

Now, by rotation we may assume that in some direction, say Ci = (1, 0, . . . , 0), 
we have by definition of the index 70 that 

70 = min{A; G N : d^F{p,Zt,u)U=o ^ 0}. 
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Then by taking A^ = 1 and y = uiin Theorem 29, we get the required estimate 
(5.23). Conditions hmsup4_^o^7o(S*) < 70 and hminf^^oo ><o(S*) > ensure 
that the dependence on t in this argument is uniform, so that the constants 
are also uniform in t. Finally, the estimate for l2{t,x) differs from that for 
Ii{t,x) in exactly the same way as in the proof of Proposition 24, so we can 
omit the repetition of the argument there. 



Suitable U'-L'^ estimates for solutions to (1) under Assumptions (Al) to 
(A4) can now be expressed as corollary of these two propositions (depend- 
ing whether all Fresnel surfaces are convex or not). 

Theorem 26 Assume (Al)i^^(^^, (A2), (A3) and (A4-)- Assume further that 
the Fresnel surfaces S|, k = 1, . . . ,m, are convex for t > t^ for some time 
to and that limsup4_^oo7(^fc) ^ 7 ™^^ hminff^oo x(Sy > 0, together with 
[{n - l)/7j < min(£i - 1,^2/2 - 2). Then solutions to (1) satisfy the LP-L'^ 
estimate 

\\U{t,-)\\Li{R") < Cpg,(l + t) ~^'i''''^'^mUo\\wP''-p{R"), 
for p E (1,2], pq = p + q, Tp = n{- — ^), and with e > arbitrary small if 



V G [0, 1) and for some e > if u = I. 



PROOF. Due to Theorem 20 solutions to (1) microlocalised to Zreg{N,u) 
have Fourier integral representations of the kind used in Proposition 24 af- 
ter multiplication with t^'^ (this multiplication makes symbol estimates uni- 
form in t as it eliminates the occurring sub-polynomial / polynomial terms 
in the estimates from Theorem 20). It remains to check the assumptions of 
Proposition 24. The phase functions (pjit,^) are averages over \j(t,C,) and 
therefore homogeneous and bounded from below and from above, in view of 
Propositions 3 and 23. Similar for the bounds of their derivatives in ^. The 
amplitudes satisfy the assumptions provided ii and £2 are big enough, i.e., 
[{n - l)/7j + 1 < min(£i,£2/2 - 1), combining (3.49) and (5.6). 

Finally, combining the microlocal estimate with Lemma 22 and Sobolev em- 
bedding theorem to estimate small times concludes the proof. 



We will give one example. If we follow Section 4.3 and consider homogeneous 
hyperbolic equations of higher order (4.13), 

D^M + Y. % „(t)DiD> = 0, BiuiO, ■) = Uj, j = 0, 1, . . . m - 1, (5.24) 

j+\a\=m 

for Qj^a £ Tui^}, satisfying the assumption of uniform strict hyperbolicity in 
combination with (4.16), then Theorem 26 applies and dispersive type esti- 
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mates depend on geometric properties of the Fresnel surfaces associated to 
the problem. Hence, if the problem is rotationally invariant all Fresnel sur- 
faces are given by spheres and we obtain 7(Sy = 2. The uniformity condition 
liminft_^oo ^(^D > is satisfied if the spheres stay bounded which is the case 
if ip{t, C,) is uniformly bounded away from zero. In this case the L^-L'^ estimate 



m— 1 m— 1 

n-l 



j2 iiiDr--'--^D^(t,-)iiL. <Cp,.(i+t)^-— j2 mr-'-wiw^-r: (5.25) 

j=0 j=0 

for 1 < p < 2,pq = p+q, r-p = n{- — -) and with arbitrarily small e > follows. 
If we drop rotational invariance, examples for 7(S^) G {2, 3, . . . 2[m/2j } can 
be constructed in analogy and give corresponding weaker decay rates. 

If some of the surfaces fail to be convex, decay rates can be much weaker. 

Theorem 27 Assume (Al)i^2, (^2)> (^3) and (A4), and assume also that 
hmsup4_>.oo7o(Sfc) < 7o with liminit-^oo x^oiX'k) > 0. Then solutions to (1) 
satisfy the W-U estimate 

\\U{t,-)\\L^(^Rr.)<Cpq^{l+t) ^o'^P i> ||f/o||w'P.'-p(R") 

for p G (1,2], pq = p + q, Tp = n{- — -), and with e > arbitrary small if 
V G [0, 1) and for some e > if u = 1. 

Remark 28 In comparison to results on scalar second order equations due to 
Reissig and co-authors, e.g., [7], [9] or [10], we observe an e-loss of decay even 
for the case z/ = 0. This is due to the C,-dependence of the term Fq in general. 
If u = the choice e = can be made in both Theorems 26 and 27 provided 
thatV^Fo{t,^) = 0. 



A The multi-dimensional van der Corput lemma 

We now give the multidimensional version of the van der Corput lemma used 
in the essential way in the proof of Proposition 24, as well as in the proof of 
Proposition 25. 

Theorem 29 ([12]) Consider the oscillatory integral 

where N > 1, and u is a parameter. Let 'y > 2 be an integer. Assume that 

(Al) there exists a sufficiently small 5 > such that x ^ C'o" (-85/2(0)), where 
-S<5/2(0) is the ball with radius 5/2 around 0; 
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(A2) $(a;, u) is a complex valued function such that Im $(a;, u) > for all x G 

supp X (^''T'd all parameters u; 
(A3) for some fixed z G supp x> the function 

F{p, u, u) := ^{z + pu, u), |w| = 1, 

satisfies the following conditions. Assume that for each fi = {u, u), function 
F{-,fi) is of class C^"^^ on suppx, and let us write its 7*'^ order Taylor 
expansion in p at as 

7 
Fip, /^) = H aj(/x)p^ + R^+iip, p) , 
i=o 

where R^+i is the remainder term. Assume that we have 
(Fl) ao(/i) = ai(/i) = for all p; 

(F2) there exists a constant C > such that I]J=2 ki(/^)l ^ C* for all p; 
(F3) for each p, \dpF{p, p)\ is increasing in p for < p < 6; 
(F4) for each A; < 7 + 1, dpF{p, p) is hounded uniformly in < p < 6 and p; 



-1, there exists a constant C^ > 



N_ 
.7 J 

such that \d^a{x, z/)| < Ca for all x G suppx and all parameters u. 



(A4) for each multi-index a of length |q;| < 



Then there exists a constant C = Cjv^-y > such that 

_N_ 

\I{^, ^)\ ^ C{1 + A) 1 for all A G [0, 00) and all parameters u. 
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